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Using both two orbital and five orbital models, we investigate the quasiparticle interference (QPI) 
patterns in the superconducting (SC) state of iron-based superconductors. We compare the results 
for nonmagnetic and magnetic impurities in sign-changed s-wave cos(A;a;)-cos(A;^) and sign- unchanged 
I cos(A;cc) • cos(A;y)| SC states. While the patterns strongly depend on the chosen band structure 
details, the sensitivity of peaks around (zbvr, 0) and (0, ±7r) wavevectors on magnetic or non-magnetic 
impurities, and on sign-changed or sign-unchanged SC orders is a common feature. Our results 
strongly suggest that the QPI can provide a direct evidence of the pairing symmetry in the SC 
states. 

PACS numbers: 74.25Jb,74.20-z 



I. INTRODUCTION 

Recently, the discovery of high temperature supercon- 
ductivity in oxypnictide compounds^i^i^i^s^ stirred great 
interests in the condensed matter community. One im- 
portant problem is to elucidate the pairing symmetry of 
the order parameter of the superconducting state. Theo- 
retically many possible gap pairing symmetries have been 
proposed. Due to the proximity of the superconducting 
state to a collinear antiferromagnetic state, a magnetism- 
based mechanism has emerged in both the weak and 
strong coupling approaches. This mechanism suggests 
that an extended s-wave pairing symmetry is favored^^^^^. 

The weak-coupling approach favors an s-wave (so 
called s±) state in which the relative sign of order pa- 
rameters changes between the hole and electron pockets. 
However, the weak-coupling approach does not specify 
the exact form of order parameter. In a recent paper^, 
we showed that, in strong-coupling, the pairing symme- 
try is determined mainly by the next nearest neighbor 
antiferromagnetic exchange coupling J2^^'^^^'^'^ and has an 
explicit Sx2y2 form in momentum space, cos{kx) ■cos{ky). 
This result is completely independent of any model, as 
long as the dominating interaction is next-nearest neigh- 
bor J2 and the Fermi surfaces are located close to the F 
and M points in the Brillouin zone. The cos{kx) ■cos{ky) 
changes sign between the electron and hole pockets in 
the Brillouin zone. In this sense, it resembles the order 
parameter, s±, proposed through general weak-coupling 
arguments^. 

The magnitudes of superconducting gaps measured by 
angle-resolved photo-emission spectroscopy (ARPES) on 
different Fermi surfaces are in good agreement with the 
simple cos{kx) • cos{ky)^^^. The magnetic properties 
in the SC state have also been shown to be consistent 
with the proposed pairing symmetr}*i^ii^iiiii^. Although 
several theoretical works ^^iSSs^iS^ propose different ways 
to measure the sign change between the electron and hole 
pockets, directly probing this change is still a fundamen- 
tal experimental challenge. Without any detailed cal- 



culations, a theoretical suggestion for probing the sign 
change through quasiparticle interference (QPI) in the 
presence of magnetic and nonmagnetic impurities, has 
been made in^. 

The QPI can be probed directly in modern STM 



experiments 



23,24 



and has been intensively studied in 



copper-based high temperature superconductors. In the 
presence of impurities, elastic scattering mixes two eigen- 
states with different momentum ki and k2 on the same 
contour of constant energy and a scattering interfer- 
ence pattern appears as a modulation in the local den- 
sity states (LDOS) at wavevector q = k2 — ki. Such 
kind of interference pattern in the wavevector space can 
be observed in the Fourier transform scanning tunnel- 
ing spectroscopy (FT-STS)^^-'^. The quasiparticle scat- 
tering between regions in the k space with high den- 
sity of states (DOS) yields peaks or arcs in the FT- 
STS. For example, in the d-wave pairing SC state, 
many QPI dispersive peaks can be identified; in the 
cuprates, they provide details of the band structure, 
the nature of superconducting gap. or other competing 
ordcrj^^^^^^^^^^^^^ 

In this paper, we perform a detailed investigation of 
the QPI in iron-based superconductors. We use both 
two orbital and five orbital models. In general, the QPI 
strongly depends on the bare band structure. The QPI 
patterns change significantly from a two orbital model 
to a five-orbital one, which suggests that the QPI can 
provide direct information of the detailed band structure 
and orbital degrees of freedom. By carefully examin- 
ing the pattern, we can also identify common features of 
the QPI pattern in both models, which are tied to the 
symmetry of SC order parameter and the impurity type. 
These general features include: (1) The intra-orbital scat- 
tering by impurities always dominates the inter-orbital 
scattering. The latter has negligible effect on the QPI 
(though it breaks discrete C4 symmetry of the patterns); 
(2) Unlike in the d-wave SC state of cuprates where a 
large density of states at the banana tips cause disper- 
sive features in the QPI^^, the nodeless s-wave has lit- 
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tie density of states inside SC gaps and hence no strong 
points dominate the scattering; (3) A magnetic impurity 
always causes a broad and large peak near q = (0, 0) in 
the QPI; this stems from intra-band scattering. For a 
non-magnetic impurity, the intensity around q = (0, 0) is 
small. This result can be used to distinguish two types of 
impurities; (4) The peaks around (±7r, 0) and (0, ±7r) are 
sensitive to both the type of impurities and to the sign 
change of the SC orders between the electron and hole 
pockets. Magnetic impurities along with sign-unchanged 
SC orders or non-magnetic impurity with sign-changed 
SC orders cause strong interference peaks. Finally, as 
in a fully-gapped s-wave SC state the results from a full 
T-matrix calculation do not differ considerably from re- 
sults of a simple first order perturbation calculation^^^^^-, 
we are able to also provide an analytic derivation of the 
above results. 



II. TWO-ORBITAL MODEL AND SINGLE 
IMPURITY SCATTERING 

We first investigate the QPI in a simple two orbital 
model^^^^'^^. The mean field Hamiltonian of the model 
in SC states is written SiS H = ^t(k)5(k)^(k) with 



B{k) = 



A^(k) 




Ai(k) 






^.(k) - 
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Nambu formalism. The single-particle bands read 



where ^'•'(k) 



ex{kx^ ky) = —2ti cos kx — 2^2 cos ky — Ats cos kx cos ky 
^y{kx •) ky) = ^x{ky J kx) J ^xy{kx •) ky) = —4:t4 sin kx sin ky 

where ti = —1, t2 = 1.3, ts = = —0.85 and fi 
is chosen in the electron-doped regime. Hereafter, \ti\ 
will be used as the energy unit. For Sx'2y2 -wave pair- 
ing, the order parameter is Ai{kx^ky) = A2{kx^ky) = 
Ao cos kx cos ky^'^^. 

The Green's function for the clean system is 

G°(k,cc;) =G^(k,k,cj) = [{u;^iS)I - B{k)]-\ (2) 

where / is the identity matrix and 5 is the energy width 
broadening. In this work, we only consider a single im- 
purity with potential ~ (5(x) so that the impurity matrix 
y(ki,k2) = F is independent of k. The impurity in- 
duced Green's function is expressed as 

5G(ki,k2,u;) = G°(ki,u;)T(ki,k2,u;)G°(k2,u;). (3) 

Standard perturbation theory gives 



V + VT^{u;)V + VT^{u;)VT^{u;)V - 

[i-vr'{u)]-'v, 
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FIG. 1: The spectral function ^(k, cj) in the unfolded Bril- 
lioun zone, for 5^,2^2 with Ao = 0.1. Darker regions corre- 
spond to larger values of A hence larger DOS in k space. 



where 



G°(k,a>). 



(5) 



Consequently, the Fourier transform of the (induced) lo- 
cal density of states is 



(4) 



27r/ 



5(k,q,a;), 



(6) 
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where q = k' — k and 

5(k,q,w) = (5Gii(k,k',w) - 5Gl^{k',k,Lo) 
SG33(k,k',Lo)-5GU'k',k,Lo). 



(7) 
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FIG. 2: Bulk density of states p for 5^,2^2 as a function of a;, 
Ao = 0.1 and no impurities. Some special values are marked 
by red dots. The energy broadening width is 5 = 0.002. 

Due to the multi-orbital nature of the band model, we 
distinguish different types of impurities. They are 

V = l^ntra = /® 0° j^y^^ (8) 

for an impurity with only intra-orbital scattering, and 

V = l^nter = ® f ^° J!.), (9) 



for inter-orbital scattering, where the upper (down) sign 
corresponds to magnetic (non-magnetic) impurity. Since 
it was argued that for cuprate superconductors this T- 
matrix method is valid when impurity scattering strength 
is much larger than the maximal pairing gap^^, we take 
Vq = 4Ao in our calculation. Our following results do 
not depend on Vq as long as Vq is much larger than Aq. 



A. Numerical results 

We first calculate electronic properties in an impurity- 
free system. Most of the results in this section have been 
already computed in less detail in^^. In Fig. [H we plot 
the spectral function 

A{k,cu) = --Im[G?i(k,c.) +G°3(k,a;)] (10) 

TT 

of the clean system at different co. The typical value of 
order parameter Aq = 0.1 will be used throughout this 
paper. It should be noted that for co and —a;, the shapes 
(topology) of the contours of constant energy (CCE) are 



almost identical, but the numerical values of A on these 
contours are remarkably different. As an example, let 
us focus on the regions near one of the M points (7r,0), 
where the CCE consists of two semi-oval-circles, or two 
complete oval-circles due to the periodicity of the Bril- 
lioun zone (BZ). As can be seen from Fig. [I](a) and (b), 
for negative cj, the spectral weight on the inner circle is 
larger than that on the outer one. The situation is op- 
posite for positive co (Fig. [T] (g) and (h)). This leads 
to different scattering interference patterns at icj, as we 
will see later on. 




FIG. 3: (color online) ^/>(q,a; = 0.07) for non-magnetic 
impurity with Mntra (black solid), Mnter (red dash) and 
Vmix = V^intra + Winter (grccu dash dot) aloug thrcc direc- 
tions: (a) (— 7r,0) (7r,0), (b) (— tt, — tt) (tt^tt) and (c) 
(-7r,7r) (tt, -tt). 

The bulk density of states p{u;) = ^(k, co) is plot- 
ted in Fig. [21 It is fully gaped within - (-0.05,0.05) 
and the coherent peak occurs at ~ ±0.076^^. 

To exemplify some scattering amplitudes, we plot 
5p{q) for a non-magnetic impurity near the edge of the 
gap {co = 0.07) along three special directions in Fig. 
[3l Two observations in these figures are common to 
all our results. First, p(q) for inter-orbital scattering 
in the two diagonal directions are quite different (com- 
pare the red dash lines in Fig. [3] (b) and (c)). This 
is not surprising since the inter-orbital scattering such 
as c[ ^ |C2,k,T breaks the symmetry between directions 
(— TT, — tt) (tt, tt) and (— 7r,7r) (tt, — tt). Second, the 
amplitude of intra-orbital impurity is stronger than that 
of the inter-orbital one; therefore the intra-orbital scat- 
tering is dominant when both are present. As such, in 
the following, we present only numerical results for intra- 
orbital impurities. 

In Fig. m and Fig. [5l we show the two-dimensional 
contour of the scattering pattern Sp{q) for non-magnetic 
and magnetic impurities, respectively. The ^p(q) profiles 
along special directions are plotted in Fig [6l In Figs. [H 
and [5l the most prominent features for all co values are 
two intersecting ovals (see also the peaks directed by red 
arrows in Fig. [6]), reflecting the strong intra-pocket scat- 
tering between equal-energy curves near M points with 
the largest DOS. These scattering processes are labeled 
by red arrows (1 and 2) in Fig. [1] (b), where scattering 
wavevectors outside the first BZ (e.g., arrow 2) should be 
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{g)uj = 0.09 (h)c^ = 0.2 

FIG. 4: Sp{q) for non-magnetic impurity with intra-orbital 
scattering, Vb = 0.4. A 200 x 200 lattice in k-space is used in 
numerical integration of equation (|5]) and the energy broad- 
ening width 6 = 0.005. 



understood as their equivalent counterparts in the first 
BZ. By increasing on the negative (positive) energy 
side, the size of the ovals decreases (increases) because 
the important scattering takes place between the inner 
(outer) circle of CCE near M with the largest DOS. At 
a definite energy uj^ this gives two peaks along direction 
M ^ r, corresponding to the scattering along the major 
and minor axis of the oval CCE respectively. Due to the 
congruence of these CCE oval circles, they always inter- 
sect on the diagonal line ( ^ F), therefore only one 
peak can be observed along this direction. For uo far out- 




(g)a; = 0.09 (h)c^ = 0.2 

FIG. 5: The same with Fig. H] but for magnetic impurity, 
Vo = 0.4. 

side the gap, there is no noticeable difference between the 
magnetic and non-magnetic impurities (compare Fig. [H 
(a) and (h) with Fig. [5] (a) and (h)). This is understand- 
able since the tendency of the (Cooper) pair breaking 
due to a magnetic impurity is most significant near the 
Fermi level. Within the single-particle scattering regime 
considered here, for \uj\ < 0.05 inside the gap, no in- 
terference pattern is expected due to lack of scattering 
states. This is confirmed (but not shown here) by the 
fact that, when decreasing the imaginary part of the en- 
ergy S in Green's functions, the peaks of (^p(q, > 0.05) 
are sharper, while (^p(q, < 0.05) vanish trivially for all 

q- 
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FIG. 6: (color online) Profiles of (5p(q,a)) along M ^ T ^ T' 
for (a) non-magnetic impurity and (b) magnetic impurity. 
The data are shifted vertically relative to each other for clar- 
ity. 



Moreover, a large peak appears around F (blue arrows in 
(Fig. m (b)). This suggests that the magnetic impurity's 
ability to localize the quasiparticle is weaker than that of 
the non-magnetic one. 
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There are additional peaks around the M points, as di- 
rected by green arrows in Fig. [6l which play an important 
role in distinguishing different types of impurities. When 
decreasing they move steadily toward M . These orig- 
inate from the inter-pocket scattering as demonstrated by 
green arrows (3 and 4) in Fig. [T] (b). The differences of 
these peaks between non-magnetic and magnetic impuri- 
ties are clear near the gap edges, suggesting strong depen- 
dence on coherent factors due to impurities, as well as on 
DOS contour of the clean system^^^^^. For non- magnetic 
impurity (Fig. [6] (a)), the peaks are much more sharper. 



FIG. 7: (color online) Same as Fig. [6] but for the case without 
sign change, Ai {kx,ky) = A2{kx,ky) = | Ao cos /c^ cos /c^ |. (a) 
Non-magnetic impurity and (b) magnetic impurity. 



B. The effect of sign change 

The Sx-2y2 pairing differs from the conventional 5-wave 
pairing that it changes sign in the BZ. To investigate the 
physical consequence of this effect, we artificially prohibit 
this sign change by letting Ai{kx,ky) = A2 (kx^ky) = 



6 



I Ao cos kx cos ky \ in Eq. ([T]) . This may not correspond to 
any realistic physical system, but can reveal, by compari- 
son, the effects of the sign change of the order parameter. 
We show the profile of Sp{q) In Fig. [3 The most observ- 
able feature is that, contrary to the sign-change case, the 
inter-pocket peak around the M point is now sharper for 
the magnetic impurity. 



C. The QPI in other pairing symmetry 

We briefiy discuss the (1^2 _y2 pairing symmetry as 
an example of a case with gapless nodal quasiparticles. 
The energy contours and interference patterns are plot- 
ted in Fig. El Fig. [9] and Fig. [lOl Besides the ro- 
bust intra-pocket scattering around F and the impurity- 
type sensitive peaks near (±7r, 0) and (0, ±7r), the most 
specific feature is the finite DOS within the pseudo-gap 
uj e (—0.05,0.05), giving rise to small but finite QPI in 
this region, as can be seen in (e), (f) and (g) of Figs. [9]and 
[TOl This originates from the intra-hole pocket scattering. 
When \uj\ 0, the QPI concentrates on the diagonal 
directions as expected from the band structures in Fig [HI 



D. Analytical analysis 

We now try to understand the above numerical results 
analytically. The two-orbital Hamiltonian can be diago- 
nalized by a unitary transformation 



C/t(k) 



where 



ex(k) 



ey(k) 



f/(k) 



ei(k) 






e2(k) 



,(11) 



U = 



cos(6lk/2) 
sin(^k/2) 



cos(6'k/2) 



and tan(6'k) 



2ex«(k) 



(12) 

The integral over BZ in calcu- 



,(k)-e^(k)- 

lating the T-matrix in Eq. ([5]) makes the analytical treat- 
ment untr act able. Fortunately we numerically verified 
that for the strength of impurity in this work {Vq = 0.4), 
the first order expansion in Eq. (|4|) is sufficiently precise 
(with error less than 2%). In the following, we safely take 
T = V. In the band-basis, the V-matrix for intra-orbital 
impurity is in the following k-dependent form 

,os(^) sin(^)\ /T/o 
.sin(^) cos(^);® V ±^0 



y(k,k') = 



as) 



where the upper (lower) sign is for the magnetic (non- 
magnetic) impurity, as in Eqs. (|8|) and (|9|). The induced 
Green's function in Eq. (j3|) can now be transformed to the 
band representation. After a lengthy but straightforward 
calculation, we obtain 
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FIG. 8: The spectral function ^(k, u) for ( 
metry with Ao =0.1. 



y2 pairing sym- 



dp{q) = [Pi(k, q) + P2(k, q)] Vo COs\ 
+ [Qi(k,q)+Q2(k,q)]Fosin2( 



2 / — ^'k-l-q 



) 

(14) 
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where Pi and Qi denote the contributions from intra- the form 
pocket and inter-pocket scattering respectively, and take 



P , (w + ei(2)(k))(a; + ei(2)(k + q))±A(k)A(k + q) 

^1(2) (k,q) = -7-2 



- A(k)2 - ei(2)(k)2)(a;2 - A(k + q)2 - ei(2)(k + q)^) ' 
^ ^ (g; + ei(2)(k))(a; + 62(1) (k + q)) ± A(k)A(k + q) 

c^i(2)l ,q; (,,2_A(k)2-ei(2)(k)2)(u;2_A(k + q)2-e2(i)(k + q)2)- ^ 



From these expressions we see that, for uo ~ ±A(/c/), 
and on the energy contour ei(k) ~ 0, the magnetic im- 
purity contribution is almost zero \i q ^ (tt, 0) for the 
sign-changing s-wave. The counterpart is true if A does 
not change sign: the contribution by magnetic impurity 
is now much larger than that by non-magnetic impurity. 
Furthermore, at ~ the contribution by magnetic im- 
purity is much larger than the non-magnetic impurity. 
These are consistent with our numerical results and pre- 
vious theoretical argument^. 

Now let us turn to the inter-orbital case. After the 
same process, one obtains 

The result is 

Mq) = [P.(k,q)+P2(k,q)]yo ^"^^''-^+;"^^''-^'^^ 
+ [Q,(k,q)+Q2(k,q)]yo ^°^^^'-^7°^^''-^'^\ (17) 

We notice that around the F point Fermi surface, 
changes from to 47r; around the M point, changes 
from to TT and then back to 0. In either case, we have 
/^^ cos{Ok)(Pk ~ 0, and therefore the inter-orbit scatter- 
ing is always much smaller than the intra-orbit scattering, 
which is also seen in Fig. [3l 



III. FIVE-ORBITAL MODEL AND NUMERICAL 
RESULTS 

A better fit to the LDA band structure in iron-based 
superconductors is given by a five orbital proposed in^^. 
To investigate the model-dependence of the scattering 
patterns, we now perform all the above calculations em- 
ploying this five-orbital model augmented by an intra- 
orbital Sx'2y2 pairing symmetry. 

In Fig. [Til we show the bulk DOS for the clean system. 
The coherence peaks appear at ±0.082 and the system is 
fully gapped within ~ (—0.05,0.05). In Fig. [121 we plot 
the spectral function v4(k, co) in the unfolded Brillioun 
zone, where the energy contours and their weight can be 
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FIG. 11: Bulk density of states p for five-orbital model with 
3^2 y2 pairing as a function of a;, Ao =0.1 and no impurities. 
The energy broadening width 6 = 0.002. 

clearly seen. In Fig. [13] and Fig. [TH we plot the in- 
terference pattern Sp{q) for non-magnetic and magnetic 
impurities, respectively. Their profiles in the direction 
M ^ F ^ F' are plotted in Fig. [151 Finally, for compar- 
ison, we also plot the results for order parameter without 
sign change in Fig. [161 

It is clear that the QPI in the five orbital model is quite 
different from that in the two orbital model. This major 
difference comes from the distribution of density of states 
at the Fermi surfaces. For example, compared with the 
QPI in the two-orbital model, scattering around F now 
has its origin in the intra-pocket scattering within the 
hole pockets - in the two orbital model it originates from 
intra-electron pocket scattering. The density of states 
is higher in the electron than in the hole pockets in the 
two band model. The opposite is true in the five orbital 
model. Another clear difference is that due to the ex- 
istence of additional orbitals, there exist square shaped 
profiles in Fig. [13] (b) and (c) which correspond to the 
scattering process labeled by arrow 1 in Fig. [12] (b); and 
the circle shaped profile in Fig. [13] (d) corresponds to 
arrow 2. These features are absent in the two orbital 
model. 

However, there are also common features in both mod- 
els. The broad and large peaks at q = (0,0) for mag- 
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FIG. 12: The spectral function v4(k, cj) in the unfolded Bril- 
lioun zone, for five-orbital model with Ao = 0.1. Darker 
regions correspond to larger values of A hence larger DOS in 
k space. 
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FIG. 13: ^p(q) for five-orbital model, non-magnetic impurity 
with intra-orbital scattering, Vb = 0.4. A 160 x 160 lattice in 
k-space is used in numerical integration of equation (|5]) and 
the energy broadening width b — 0.001. 



netic impurity appear in both models. More importantly, 
the (±7r, 0)/(0, ±7r) sensitiveness on magnetic or non- 
magnetic impurity, and sign change remains the same. 
For example, when A changes sign, the peak around 
q = (±7r, 0)/(0, ±7r) for non-magnetic impurity (Fig. [15] 
(a)) disappears in the case of a magnetic impurity (Fig. 
[T5l (b)). On the contrary, when A does NOT change 
sign (Fig. [T6|) . the peak at (±7r, 0)/(0, ±7r) is related to 
magnetic impurity. The sensitiveness of the interference 
pattern around (±7r, 0)/(0, ±7r) corresponds to the inter- 
pocket scattering labeled by arrows 3 and 4 in Fig. [12] 



(b) and has been explained explicitly in the two-orbital 
model. Eqs. ([T4|) and ([T5|) and the arguments follow- 
ing them do not depend on the number of bands and 
therefore these features are rather universal. Moreover, 
the (±7r, 0)/(0, ±7r) sensitiveness on the order parameter 
sign change (Fig. [T6|) is quite similar with the two-orbital 
model. 
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(c)cj : 



-0.08 



(e)cj = 0.07 




{h)u; = -0.09 




-0.07 



(f)u; = 0.08 



I + 



(g)c^ = 0.09 (h)cj = 0.2 

FIG. 14: Same as Fig. [T3]but for magnetic impurity, Vo = 0.4. 

IV. CONCLUSION 

In summary, we have investigated in detail the struc- 
tures of the QPI in iron-based superconductors within 
the current available two orbital and five orbital models. 
The results obtained here suggest that the QPI can be 
used to determine the band structure and orbital degrees 
of freedom in these materials and can also provide evi- 
dence of the SC pairing symmetries. In this calculation, 
we have ignored possible three dimensional effects, more 
relevant for the 122 materials^i^. The physics associ- 




ated with the third dimension and possible competing 
orders or coexistence states will be addressed in the fu- 
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FIG. 15: (color online) Profiles of 5y9(q, u) along M ^ F ^ F' 
for (a) non-magnetic impurity and (b) magnetic impurity in 
five-orbital model. 
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